Wave function correlations on the ballistic scale: 
Exploring quantum chaos by quantum disorder 
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We study the statistics of wave functions in a ballistic chaotic system. The statistical ensemble 
is generated by adding weak smooth disorder. The conjecture of Gaussian fluctuations of wave 
functions put forward by Berry and generalized by Hortikar and Srednicki is proven to hold on 
sufficiently short distances, while it is found to be strongly violated on larger scales. This also 
resolves the conflict between the above conjecture and the wave function normalization. The method 
is further used to study ballistic correlations of wave functions in a random magnetic field. 
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Introduction. Understanding of statistical properties 
of eigenfunctions of a quantum system whose classical 
counterpart is chaotic and their relation to the underly- 
ing classical dynamics is one of the key problems studied 
in the field of quantum chaos. Among various applica- 
tions, wave function correlations are important for statis- 
tics of electron transport through quantum dots, see [ffl 2 ] 
and references therein. It was conjectured by Berry 3] 
that an eigenfunction of a classically chaotic system ( "bil- 
liard" ) can be represented as a random superposition of 
plane waves with fixed absolute value k of the wave vector 
(determined by the energy k 2 /2m = E, where m is the 
mass and we set K = 1). This implies Gaussian statistics 
of the eigenfunction amplitude 
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determined solely by the correlation function 
sider a two-dimensional system) 



C(r 1 ,r 2 ) = (^*(r 1 )V(r 2 )) = J (*|n - r a |)/V. (2) 

Here (3 — 1 or 2 for a system with preserved (respec- 
tively broken) time reversal symmetry, V is the system 
area, and Jo the Bessel function. For defmiteness, we will 
consider the case (3 = 2 below; generalization to systems 
with P = 1 is straightforward. 

Hortikar and Srednicki ||| noticed that Eqs. (|l|), (||), 
which do not depend on any details of the dynamics, 
may only be valid for sufficiently small spatial separa- 
tion. They generalized Berry's hypothesis and conjec- 
tured the Gaussian statistics ([!]) with a more general, 
system-specific kernel C(ri,r2) replacing Eq. (0), 
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where G sc is a semiclassical Green's function Q. This 
proposal was supported by the observation JBj that the 
result obtained in Ref. Q for two-point correlations in a 
diffusive system is consistent with the Gaussian statistics. 



The conjecture (||), (g), while physically very appeal- 
ing, obviously requires a formal derivation. Furthermore, 
when taken literally, this conjecture contradicts the wave 
function normalization, 



dr [(^W(r')l} - (|^ 2 (r)l><I^V)l>] = 0, (4) 



since the integrand is equal to C 2 (r, r') > according to 
(0). Therefore, limits of validity of this conjecture have 
to be understood. All this points to a need in a system- 
atic study of wave function statistics in ballistic systems, 
which is the aim of the present paper. 

Eigenfunction statistics in a ballistic system. In or- 
der to speak about the wave function statistics V{ip(r)}, 
one should first define an ensemble over which the av- 
eraging goes. Such an ensemble can be generated @] 
by adding to a system under consideration a random 
potential U(r) characterized by a correlation function 
W(r — r') = (U(r)U(r r )) with a correlation length d. 
Parameters of this random potential are assumed to sat- 
isfy A; -1 <C d <C l s *C L <C hr, where l s (lt r ) is the 
single-particle (respectively transport) mean free path, 
and L is the characteristic size of the system. The con- 
dition l tr ^> L ensures that the additional disorder does 
not influence the classical dynamics of the system, while 
the inequality l s <C L guarantees that the ensemble of 
quantum systems is large enough to produce meaningful 
result. Note that the potential is smooth, kd ^S> 1, since 
hr/l s ~ (kd) 2 . On the technical side, introducing the 
additional random potential allows us to apply, with a 
proper generalization, methods developed earlier for dif- 
fusive systems (see M for a review). 

After the ensemble averaging, the problem is described 
by a ballistic a- model of a supermatrix field Q(r, n) with 
the action 
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where H is the Hamiltonian of the system (without dis- 
order), v is the density of states, w(n, n') = 2-Ki>W(k\n — 
n'|) is the scattering cross-section by the random poten- 
tial, and n is a unit vector characterizing the direction of 
velocity on the energy surface. On the scales 3> l s Eq. (||) 
reduces to the form proposed in Refs. Jl2] , ^3| . The two- 
point correlation function of the wave function intensities 
is expressed in this approach as 



(|V 2 (ri)^ 2 (r 2 )|> = Urn ^([G 1 i(ri,r 1 )G 22 (r 2 ,r 2 ) 

17— >Q 7T 

+ Gi 2 (ri,r 2 )G 2 i(r2,ri)]) 5[ Q], (6) 

where A is the mean level spacing, G is the Green's func- 
tion in the field Q, 
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and the subscripts 1,2 refer to the advanced-retarded 
decomposition (the boson-boson components being im- 
plied). We first evaluate Eq. @ in the zero- mode ap- 
proximation, Q(r) = Qo- The Green's function ([?]) is 
given in the leading order by 

G (r, r') = ilmG R (r, r')Qo + ReG R {r, r'), (8) 
G R (r, r') = {r\(E-H + i/2r s )- l \v') . (9) 

If the points r, r' are separated by a distance ^> l s from 
the billiard boundary, Eq. (^|) reduces to G R (r, r') = 
J (fc|r - r'|)e-l r - r 'l/ 2Zs . Substituting Eq. (|) in (§) and 
expanding the action (||) up to the linear-in-77 term, 
S[Q] ~ irvrjVStrQoA, one finds, in full analogy with the 
case of diffusive systems, 

^(^(^(raJD-l + fc,^,^); (10) 
k q {v, r') = ImG fl (r, r')ImG fl (r', r)/(7ri/) 2 , (11) 

with the two contributions on the r.h.s. of ( |i~0| ) originat- 
ing from the terms (GnG 22 ) and (Gi 2 G 2 i) in @, re- 
spectively. The result ( |l0| ) , corresponding exactly to the 
conjecture (Q), (0) of the Gaussian statistics, is in con- 
flict with the wave function normalization, as explained 
above. 

To resolve this problem, we evaluate the term (GnG 22 ) 
more accurately by expanding the Green's function (|7]) 
to the order 77 and the action (||) to the order i] 2 . While 
these terms (usually neglected in the er-model calcula- 
tions) are of the next order in 77 and may be naively 
thought to vanish in the limit 77 — > performed in @, 
this is not so, since Qq oc As a result, we get in the 
zero-mode approximation 



^ 2 (|V 2 (ri)V 2 (r 2 )|)zM-l 

= fc g (rx,r 2 ) - fcg(ri) - k q (r 2 ) + k q 



(12) 



(terms of still higher orders in r\ produce corrections small 
in the parameter Ar s <C 1), where 



k q (r)=V- 1 J d 2 v'k q {r,v'), 
k q = V- 2 [ d 2 rd 2 r'k q {r,r'). 



(13) 

The contribution of non-zero modes is found to be |?J 

^ 2 (|V' 2 (ri)V' 2 (r2)|)NZM =n B (ri,r 2 ), (14) 

where iTs(ri,r 2 ) = IIs(ri,r 2 ) — IT^(ri,r 2 ) describes 
the (integrated over direction of velocity) probability of 
classical propagation from ri to r 2 , 

ns(ri,r 2 ) = J J dnidn 2 X>(rini,r 2 n 2 ), 

CV = (tit/) -1 [8(ti ~ r 2 )S(n 1 - n 2 ) - V' 1 ] , (15) 

with the contribution II^(ri,r 2 ) of direct propagation 
(before the first event of disorder scattering) excluded. 
The symbol £ in (15) denotes the Liouville operator char- 
acterizing the classical motion Jl4| . 

We analyze now the total result given by the sum 
of ([l2|) and (|l4|). First of all, we stress that it satis- 
fies exactly the constraint (^) of wave function normal- 
ization. Next, we consider sufficiently short distances, 
|ri — T2| -C Z s - In this case the correlation function is 
dominated by the first term in the r.h.s. of Eq. jl^), 
returning us to the result (10). Furthermore, we can 
generalize this result to higher correlation functions, 



(r (riMrO • . . r(r n )Hr' n )) = - 2V{n _ iy 



S[Q] 



where the summation goes over all permutations a of the 
set {1, 2, . . . , n}, pi = 1 for % = 1, . . . , n — 1, and p„ = 2. 
If all the points , are within a distance <C l s from 
each other, the leading contribution to this correlation 
function is given by the zero-mode approximation with 
higher-order terms in rj neglected [i.e. by the same ap- 
proximation which leads to Eq. (fffl)], yielding 



V"<^(ri)V(ri) ■ ■ ■ r{r n )Wn)) = E II M^r^), 
/ F (r,r') = -ImG fl (r,r')/(vr i ,). (16) 

This result is identical to the statement (Q) of the Gaus- 
sian statistics of eigenfunctions. We have thus proven 
that the conjecture of Refs. ||, Q holds within a spatial 
region of an extension <C l s , with the kernel G(r l5 r 2 ) 
given by Eqs. @, ®. 

We turn now to the behavior of the correlator 
(|V' 2 (ri)V' 2 (i' 2 )|) at larger separations |ri — r 2 | ^> l s - In 
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this situation, the correlations are dominated by the con- 
tribution ( |l4| ) of non-zero modes. Let us further note that 
the smooth part of the zero-mode contribution ( |l2] ) (i.e. 
with Friedel-type oscillations on the scale of the wave 
length A = 2-rt/k suppressed) is exactly equal to 11^. 
Therefore, the smoothed correlation function is given by 
the classical propagator, 



y 2 (|V' 2 (ri)V' 2 (r2)|) sm ooth - 1 = n B (n,r a ), 
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independent of the relation between |ri — r2 1 and l s . The 
mean free path l s manifests itself only in setting the 
scale on which the oscillatory part of (|-0 2 (ri)^ 2 (r2)|) 
gets damped. 

A result for the variance of matrix elements related to 
Eq. ( |l7j ) was obtained in |t5| by a semiclassical method. 
Note, however, that the semiclassical treatment of [[l5| is 
only justified if one introduces a sufficiently large level 
broadening 77 ^ A, while calculating statistical proper- 
ties of a single eigenf unction requires the limit r\ <C A, 
see Eq. (||). Also, the result of |15| does not contain the 
contribution corresponding to the term V^ 1 in (fi"5|), so 
that it violates the wave function normalization in the 
same way as the conjecture (|l|), (||). 

Since we have shown that for l s <^ L the applica- 
bility of the Gaussian statistics (Q), (|J) is restricted to 
the scales <gC l s , one may be tempted to ask whether in- 
creasing l s beyond L would be favorable from this point 
of view. The answer is no; in contrast, for l s > L a 
further increase of l s reduces the range of applicability 
of the Gaussian statistics. Indeed, it is not difficult to 
show that for l s 3> L the Green's function (||) has the 
form G fl (ri,r2) ~ Jo(fc| r i — r 2|) (we assume for sim- 
plicity that the points ri , r2 are sufficiently far from the 
boundary) only for |ri — r 2 | <C l s , where l s — L 2 /l s . 
At larger distances, |ri — r 2 | > / s , the Green's function 
shows irregular oscillations with a characteristic ampli- 
tude |G fl (ri,r 2 )| ~ (kls)^ 1 / 2 independent of |i"i — r 2 |, 
which are physically due to the interference of waves 
multiply reflected from the boundary. Therefore, only 
at |ri — r2 1 <C l s the first term in ( |l2| ) will dominate and 
the statistics will be Gaussian. 

Random magnetic field. In the above we studied the 
wave function statistics of a given chaotic system by gen- 
crating an ensemble of quantum systems with the help 
of an additional random potential. Now we use the same 
approach to study the wave function statistics in a ran- 
dom magnetic field (RMF). In this case, an ensemble is 
defined from the very beginning and the introduction of 
additional weak disorder may be considered as a technical 
trick, the reason for which is explained below. 

We consider a white-noise RMF B(r) with the corre- 
lation function 



(B(r)B(r')) = TS^ (r - r'), T < fc 2 



(18) 



free path l tr = 4fc/r this problem is described by the 
conventional unitary-class diffusive cr-modcl [ ]ic| so that 
the results obtained for diffusive systems Q apply. We 
will be interested, however, in wave function correla- 
tions on much shorter - ballistic - length scales. Specif- 
ically, we will study how the Friedel-type oscillations in 
( IV^ 2 ( r i )"0 2 ( r 2 ) I } decay with increasing |ri — r 2 |. [The 
smooth part is simply V 2 (\ip 2 (ri)ip 2 (r2)\) = (irk Y\ — 
r 2 |) _1 for all |ri — r 2 | <C hr, L, as follows from (|l7)]. 

In the case of a random potential the scale for the van- 
ishing of oscillations is set by the single-particle mean 
free path l s . An attempt to get an analog of this result 
by deriving directly the ballistic a-model via averaging 
over the RMF fails, since the equation for l s obtained 
within the self-consistent Born approximation (SCBA) 
leads to an infrared-divergent and gauge-dependent re- 
sult |l6). This is a manifestation of the fact that in the 
case of a RMF the single-particle relaxation rate depends 
on geometry of the problem. 

To overcome this problem, we add an additional weak 
random potential with the mean free path much 
longer than the length scale of interest set by the RMF 
(which we will find below). Averaging over this ran- 
dom potential, we derive the cr-model in a given real- 
ization of the RMF. As explained above, the two-point 
correlation function of eigenfunction intensities on a scale 
|ri - r 2 | < lf F is given by Eqs. ([(J, (|ll|). Therefore, 
the desired oscillatory contribution reads 

(fc° sc (r 1 ,r 2 )) R MF 

= -(7n,)- 2 Re (G R (r u r 2 )G> 2 , ri)) RMF , (19) 

where G R — (E — H + i/2rJ lP ) _1 is the Green's function 
in a given realization of the RMF, and (. . .)r,mf denotes 
averaging over the RMF realizations. This (second) av- 
eraging can be performed with use of the path integral 
formalism JlJ] . The product of the two Green's functions 
in ([L9|) can be written as 



poo <.ri(Ti)=R 

(G fl (R,0)G R (0,R)) RM F= / dT x dT 2 / Vr x Vv 2 

J0 J ri (0)=0 

x exp[i(E + i/2r s RP )(T 1 + T 2 ) + iS kin - S RMF ], (20) 



where 5km — Jq Ti dt mr 2 /2 + J Q 12 dt mr|/2, and .Srmf = 
rs no /2, with s no denoting the non-oriented area between 
the two trajectories ri(t) and V2(t). The integral ( pp| ) is 
dominated by the pairs of paths being close to each other 
and corresponding to an almost uniform and straight 
motion from to R. To make this explicit, it is use- 
ful to perform the change of variables |17|, introducing 
r = n r 2 , p = (n + r 2 )/2, t+ = (T a + T 2 )/2, and 
t- = T 2 - Ti (t + > t_). The RMF-induced part of the 
action takes then the form 



r T 2 



and assume that the size of the system, L, is sufficiently 
large. On length scales longer than the transport mean 
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where v = k/m is the particle velocity, and we have split 
r into components parallel (r") and perpendicular (r ± ) 
to p w R/t_|_. While the integrals over p and r'> are es- 
sentially the same as for a free particle, that over r has 
the form of the Feynman integral for a one-dimensional 
particle in the potential z(t>r/2)|a;|. The corresponding 
Green's function g(x,x',t) reads in the frequency repre- 
sentation at x = x 1 = (which is what we need in view 
of the boundary conditions on r ) , 
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where Ai(z) is the Airy function, and 
r = (4m/r\ 2 ) 1 /3. 



(22) 



(23) 



This leads to the following result for the oscillatory part 
of the wave function correlation function 



y 2 (^ 2 (ri)^ 2 (r 2 ) 
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/RMF 



X < 



1 

■nkr 
sin(2fcr), 

sin(2fer + \Q Q \r/2l Q - vr/12) 
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— 1.05 is the lowest zero of 



where r = |ri — r 2 | 
Ai'(z), and ?o = vtq. 

We thus find that the oscillations are suppressed 
on the scale ~ l = vt = (4fc/r 2 ) 1 / 3 . Note that 
Iq is parametrically different from both the transport 
mean free path l tr = 4/c/F and the length IdHvA = 
vTdHvA = (27r/r) 1 / 2 characterizing damping of de Haas - 
van Alphen magnetooscillations of the density of states, 
Pose oc exp[— (7r/w c TdHvA) 2 ] | IT | . Difference between Zo 
and ZdHvA (in the case of a random potential both these 
scales are set by l s ) illustrates the already mentioned de- 
pendence of the single-particle relaxation rate on the ge- 
ometry of the problem in the case of RMF. 

The scale Iq was obtained in |l8| from consideration of 
certain Green's function with an obscure physical mean- 
ing. We have demonstrated that the length l determines 
an observable quantity - the scale of decay of the oscil- 
latory part of the wave function correlation function. 

Conclusions. We have studied the wave function statis- 
tics in a chaotic ballistic system. The corresponding sta- 
tistical ensemble is defined by adding a smooth random 
potential, satisfying l tr L ^> l s . The first inequality 
guarantees that the random potential does not change 
the classical dynamics, while the second one ensures that 
the ensemble of quantum systems is sufficiently large and 
provides meaningful statistics. By using the ballistic a- 
model approach we have proven that the conjecture of 



Gaussian fluctuations of wave functions [|3|J4[] holds on suf- 
ficiently short distances |r,* — rj| <C l s , while it is strongly 
violated on larger scales. Our results solve, in particular, 
the problem of inconsistency of the conjecture of Gaus- 
sian statistics with the wave function normalization. 

We have further applied these results to study the de- 
cay of Friedel-type oscillations in the correlation function 
(|^ 2 (ri)-0 2 (1*2) |) in a RMF. In this case averaging over an 
additional weak random potential yields Gaussian fluc- 
tuations of wave functions in a given realization of the 
RMF. The remaining averaging over the RMF realiza- 
tions performed via the path integral formalism leads to 
the result (pij). The scale lo for the decay of oscillations 
(playing the role of the single-particle mean free path l s ) 
is given by Eq. (E3I), providing physical meaning to a 
length found in JL8| from some formal consideration. 
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